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Abstract 

A non-perturbative expansion method which gives a well-defined analytic continuation 
of the running coupling constant from the spacelike to the timelike region is applied to 
the inclusive semileptonic decay of the r-lepton. The method allows us to evaluate R T 
by integration over the non-perturbative physical region, rather than by using Cauchy's 
theorem, and hence to incorporate threshold effects in a very direct way. Within our 
framework the difference between the effective coupling constants in the timelike and 
spacelike domains can be substantial and is not simply a matter of the standard tt 2 terms. 



The t decay process with hadronic final states represents an important test of quantum 
chromodynamics. Due to the inclusive character of the process, the ratio 

r( t~ -»• u T hadrons(7) ) 
T(r ^z/ T e z/ e (7j) 

where (7) denotes possible additional photons or lepton pairs, is a very convenient quantity 
both for a theoretical investigation and for the definition of the QCD coupling constant a s (M 2 ). 
A detailed theoretical analysis of this problem has been given in Ref. [1] ( see also Refs. [2-5], 
in which different aspects of the problem are discussed ). 

The starting point of the theoretical analysis is the expression 



where 



N 

= 77- [ n (s + ie) - n(s - ie) 

Z7T1 



n(s) = E |v^f(iw( s ) + n ug , A (s)) . (3) 



q=d,s 



The normalization factor N is defined so that in zeroth order perturbation theory R^} rt = 3. 
In the framework of standard perturbation theory the integral (2) cannot be evaluated directly 
since the integration region in (2) includes small values of momentum for which perturbation 
theory is invalid^ Instead of Eq. (2), the expression for R T may be rewritten, using Cauchy's 
theorem, as a contour integral in the complex s-plane with the contour running clockwise around 
the circle \s\ — M%. It seems that this trick allows one to avoid the problem of calculating the 
nonperturbative contribution, which is needed if one uses Eq. (2). However, the application of 
Cauchy's theorem is based on specific analytic properties of IT(s) or the Adler D function 

D(q 2 ) = <? 2 (-^) Nn 0? 2 )- (4) 

The function D(q 2 ) is an analytic function in the complex g 2 -plane with a cut along the pos- 
itive real axis. It is clear that the approximation of the D-function by perturbation theory 
1 In Ref. [6], the integral (2) has been calculated within the method of optimized perturbative series [7]. 



breaks these analytic properties. For example, the one-loop approximation for the QCD run- 
ning coupling constant has a singularity at Q 2 = Aq CD , the existence of which prevents the 
application of Cauchy's theorem. Moreover, to define the running coupling constant in the 
timelike domain, one usually uses the dispersion relation for the D function derived on the 
basis of the above-mentioned analytic properties. In the framework of perturbation theory, this 
method gives the so-called 7r 2 -term contribution which plays an important role in the analysis 
of various processes [8-12]. However, the same problem arises: the perturbative approximation 
breaks the analytic properties of A eff (g 2 ) which are required to write the dispersion relation. In 
addition, there is the problem of taking account of threshold effects. As follows from Eq. (2), 
the initial expression for R T "knows" about the threshold; for example, the value R T written in 
the form of Eq. (2) contains the information that the CKM mixing matrix element V us = for 
s < (m u + m s ) 2 . But all the threshold information is lost if one rewrites Eq. (2) as a contour 
integral and uses a fixed number of flavours for the calculation of R(s) on this contour. 

In this note we will concentrate on both aspects of the problem. Our considerations are 
based on a method of construction of the effective running coupling constant which respects 
the above-mentioned analytic properties [13,14]. In Ref. [15] we have demonstrated that in 
the framework of this approach there exists a well-defined procedure for defining the running 
coupling in the timelike domain which does not conflict with the dispersion relation. We will 
use the following definitions: A eff = aQCD/(4vr) is the initial effective coupling constant in 
the t-channel ( spacelike region ) and Xf 1 is the effective coupling constant in the s-channel 
( timelike region ). From the dispersion relation for the D-function we obtain [15] 

A e V) = -? 2 f-^Af( S ). (5) 
Jo [s — q l ) 

Thus, the initial running coupling constant \ cS (q 2 ) is an analytic function in the complex 
g 2 -plane with a cut along the positive real axis. This function does not exist for real positive q 2 , 
so the definition of the running coupling constant in the timelike domain is a crucial question. 
Here we use the standard definition of A^ ff (s) in the s-channel based on the dispersion relation 
for the Adler D-function. In this case, parametrization of timelike quantities, for example 
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R e + e -(s) or R(s), by the function A^ ff (s) is similar to parametrization of spacelike processes by 
the function A cff (g 2 ). 

The inverse relation of Eq. (5), given the analytic properties of A cff (g 2 ), is of the form 



K S (s) = -^- ^-A e V), (6) 



JL r" 'hf \eff,„2 

27ri Js-ie q 2 

where the contour goes from the point q 2 = s — ie to the point q 2 = s + ie and lies in the region 
where X eS (q 2 ) is an analytic function of q 2 . Equation (6) defines the running coupling constant 
in the timelike region which we must use to calculate R(s) in Eq. (2). As follows from Eqs. (5) 
and (6), there is a connection between the asymptotic behaviours of \ eS (q 2 ) and A^ ff (s). If the 
function A cff (g 2 ) has the asymptotic behaviour 

b \n(-q 2 /A 2 ) 

the function A^ ff (s) behaves like 

Af (S) ~ b Ms/ A 2 ) aS 56 + °° 

and vice versa. 

To write Eq. (6), it was important that the function A cff (g 2 ) had the above-mentioned 
analytic properties. For example, to use the one-loop approximation, one needs to modify its 
infrared behaviour at Q 2 = A 2 in an ad hoc manner so that the singularity at Q 2 = A 2 is absent 
in the new expression for \(Q 2 ) . A self-consistent formulation of the analytic continuation 
problem is, however, possible within the scope of a systematic non-perturbative approach. 

Let us rewrite Eq. (2) in the form 

R T = 2 f 1 dx(l - xf{l + 2x)R(M*x) , (7) 
Jo 

where R(M 2 x) can be expanded as the series 

R{M 2 T x) = r [l + n A S (M 2 x) + r 2 A 2 (M 2 x) + ■ ■ •] 

= r [1 + n\f(M 2 x)]. (8) 
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Here we have introduced the effective coupling constant in the s-channel, which incorporates 
higher-order corrections. As has been mentioned above, any approximation of the quantity R 
by a finite number of terms of the perturbative series does not allow us to calculate the integral 
(7) directly because this expression includes the integration region of small s where perturbation 
theory is invalid. In Ref. [15], the s-channel effective coupling constant has been constructed 
within the nonperturbative method proposed^ in Refs. [13,14]. This approach to QCD allows 
one to systematically investigate the low energy structure of the theory and evaluate the integral 
(7) directly without the use of Cauchy's theorem. The method is based on a new expansion 
parameter a which is connected with the original coupling constant by the following equation 



a 2 



where C is a positive parameter which we will fix on the basis of meson spectroscopy [14] . From 
Eq. (8) it is easy to see that the expansion parameter a is bounded by < a < 1 for all values 
of the initial coupling constant in < A < oo. 

The Q 2 -evolution of the parameter a is defined by 



M = /(oo) + ^ |j , (10) 

where a is the value of the parameter a at some normalization point Qo and to order a 3 the 
function f(a) has the form 

2 6 „o, 18 1 624 i h \ 5184 , H 9 , 

/ « = 481na- — + — In 1 - a + — - In 1 + - a . 11 

ar a 11 1 — a 121 121 2 

The effective coupling constant in the s-channel can be written as 

K{s) = ^Im0(a + ) , (12) 



where 



72 1 318 , , . 256, 9. 

m = -4 In a - -— + - 1„(1 - .) + - ln(l + -a) , (13) 



In Ref. [16] the method has been applied to describe the process of e + e annihilation at low energies. 
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and the value of a + obeys the following equation: 



f(a + ) =/( ao ) + ^[ln^ + i7r]. (14) 



In the calculation of the integral (7) the number of active quarks is different in various 
regions of the integration. One usually applies the matching procedure in the Euclidean t- 
channel, changing the number of active quarks at the threshold Q = £ m q with a matching 
parameter 1 < £ < 2 [17] ( see also Refs. [18,19] ). However, within the framework of a 
massless renormalization scheme, a more natural way is to introduce the threshold matching 
in the physical s-channel, where, at least in leading order, the change in the number of active 
quarks is an obvious fact associated with the energy threshold of the quark pair production. 
Moreover, any matching procedure for the coupling constant in the t-channel, for which one 
uses a condition of the type ReQ 2 > £ 2 m 2 , leads to a violation of the analytic properties of 
A eff (g 2 ), whereas the matching procedure in the s-channel maintains the analytic properties of 
the effective coupling constant X eS (q 2 ) in the complex g 2 -plane. 

By using the method under consideration we can ensure that both A s (s) and its derivative 
Ag(s) are continuous at various threshold points Sj. The set of corresponding equations for the 
parameters and Oq is as follows: 

' Im^flf 11 ) = -i-Imrf) , (15) 



bo(f-l) ^ + ' b (f) 

^Im[(at 1} ) 2 (1 + 3ar i} )] = Im [ {aff (1 + 3a ifh 



For the number of active flavours / = 3 we use the value = 4.1 obtained from the 
phenomenology of meson spectroscopy [14]. Then for / 3 we find the parameter using 
Eqs. (15). We use for the matching parameter £ the value £ = 2, but the value of R T from 
Eq. (7) is practically independent^ of the parameter £ in the interval 1 < £ < 2. 

In the alternative approach to the evaluation of R T using Cauchy's theorem one proceeds 
first by an integration by parts to convert R into D, then represents the discontinuity as a 
3 The corresponding uncertainty is smaller than 0.01% . 
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contour integral and finally opens up the contour to the unit circle in the z plane. In this way 
R T is expressed in terms of A eff as 

*r = T^i , - (1 - z f (! + *) D(M 2 z) , (16) 
2iri J\z\=i z 

where the D function has the following form 

D(M 2 z) = d [ 1 + d x \ cS (M 2 z) ] , (17) 

and to order 0(a 3 ) [] 

\ cS (Q 2 ) = ia 2 (l + 3a), (18) 
where a = a{Q 2 ) is found from Eq. (10). We can rewrite Eq. (16) in the form 

R T = i?[ 0) (l+ AR T ) , (19) 

with i4 0) defined as 

4 0) = 3(\V ud \ 2 + |K,| 2 )^ew, (20) 

where the electroweak factor Sew an d the CKM matrix elements \V u d\, \V US \ are Sew 
1.1)1 1)1 . \V ud \ = 0.9753, \V US \= 0.221 taken from Ref. [1]. 
For AR T one finds 

ARr = d x / — (1 - z) 3 (1 + z) A eff (M r 2 z) . (21) 

27T1 ^|«|=1 Z 

By using the parametrization z = — M 2 e %e and with d\ = r\ = 4 we get 

AR T = - r d9(l + 2e w - 2e 3W - e 4id ) X eS (M 2 e id ) . (22) 

7T J-iT 

It should be stressed that we can use Eq. (22) only for a fixed number of quark flavours. 

Indeed, we have applied the physical matching procedure in the s-channel and have obtained 

an analytic function A eff (q 2 ) with a cut along the positive real axis. This function "knows" , 

in principle, about all quark thresholds due to Eq. (5). For a function A cff (q 2 ) obtained in 
4 In 0(a 2 ) we obtain similar results, in accordance with the mechanism of induced convergence (for a more 

detailed account see Ref. [15]). 



such a way it is impossible to split the region of Q 2 = — q 2 into subregions in which the 
number of active quarks is fixed. The situation is similar to the situation that arises in the 
mass-dependent momentum renormalization scheme. In this scheme, all quarks contribute to 
the effective coupling constant; however, for Q 2 « M 2 the quark with mass M q is irrelevant 
because its contribution to the running coupling is very small. 

To check the consistency of our method let us first consider a fixed number of active quarks 
/ = 3 and use R T = 3.56 as an input. In this case, using Eqs. (7),(8),(12)-(14) we can find the 
parameter ao and then verify that Eqs. (7) and (16) give the same results for R T : R T = 3.560 
as should be the case. 

To take into account the threshold effects, we have used Eqs. (15) to find the parameter 
(jU) f or f -L 3 anc i the conditions that the CKM matrix elements V u d = for s < (m u + rrid) 2 
and V us = for s < (m u + m s ) 2 . We use here the following values of the quark masses: 
m u = 5.6 MeV, = 9.9 MeV and m s = 199 MeV. In this way, we obtain from Eq. (7) for R T 
the value 3.552 instead of 3.560. One can see that the threshold effects for R T are about 0.2 %. 

Now, taking the experimental value R T = 3.552 [20] as an input, we obtain a s (M 2 ) = 0.37 
and a(M 2 ) = 0.40. The values of the coupling constant in the s- and t-channels are clearly 
different from each other; the ratio is a s (M 2 )/a(M 2 ) = 0.92. 

The experimentally measurable quantity R T can be parametrized both by the function a s (s) 
defined in the time-like region and entering into the initial expression for R T (Eq. (2)) and by 
the running coupling constant a(q 2 ) used in the contour integral (Eq. (16)). The perturbative 
expansion does not allow one to perform the integration in Eq.(2) directly because it involves 
a non-perturbative region. Instead, one usually uses the perturbative formula to evaluate the 
contour integral (16). However, we believe this to be inconsistent because the analytic properties 
which are required to write down the Cauchy integral are not respected by the perturbative 
formula. 

In the present paper we have proposed a method which allows one to evaluate both the 
initial integral for R T and the expression obtained by the use of Cauchy's theorem. They are 
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equal in the absence of threshold effects, which, however, can be easily incorporated in the first 
method. We have also demonstrated that the distinction between the functions a s (s) and a(q 2 ) 
is not simply a matter of the standard n 2 terms, which may be important for understanding 
certain discrepancies [19, 21] arising in the determination of the QCD coupling constant from 
various experiments. For a more detailed quantitative analysis of the experimental data using 
our method it will be necessary to include higher-order corrections. In the near future we plan 
to extend the calculation to 0(a 5 ). 

The authors would like to thank D.I. Kazakov and A.N. Sissakian for interest in the work 
and useful comments. 
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